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AN EXTENSION OF MACMAHON’S EQUIDISTRIBUTION 
THEOREM TO ORDERED MULTISET PARTITIONS 


ANDREW TIMOTHY WILSON 


Abstract. A classical result of MacMahon states that inversion number and 
major index have the same distribution over permutations of a given multiset. 
In this work we prove a strengthening of this theorem originally conjectured 
by Haglund. Our result can be seen as an equidistribution theorem over the 
ordered partitions of a multiset into sets, which we call ordered multiset par¬ 
titions. Our proof is bijective and involves a new generalization of Carlitz’s 
insertion method. This generalization leads to a new extension of Macdonald 
polynomials for hook shapes. We use our main theorem to show that these 
polynomials are symmetric and we give their Schur expansion. 
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1. Introduction 

1.1. Permutation statistics. Given a composition a of length n (i.e. a vector 
of positive integers of length n), we let ©^ be the set of all permutations of the 
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multiset {i“* : 1 < z < n}. For a permutation a G ©a, written in one-line notation, 
the descent and ascent sets of cr are 

Des(cr) = {i : (Tj > a^+i} Asc{a) = {i : ai < ai+i} 

The inversions of a are the pairs 

Inv((T) = {(z, j) ■■ 1 < i < j < n, cr* > aj}. 

It will be convenient to refine the set of inversions in the following manner: 

Inv*’° = {{i,j) - i < j <n, cr* > aj} Inv°’^ = {{i,j) : 1 < z < j, cn > (Jj}. 

These are the elements of Inv((T) whose first (resp. second) coordinate is z (resp. j). 
These sets allow us to define several statistics on ©„: 

des((T) = I Des((T)| asc((j) = | Asc((t)| inv((T) = | Inv((T)| maj(cr) = ^ z. 

zGDes((7) 

These statistics are known as the descent number, ascent number, inversion number, 
and major index of a, respectively. We will also make use of two refinements of 
inversion number: 

inv*’'^(cr) = I Inv*’'^(cr)| inv'^’-^(cr) = | Inv'^’'^(cr)|. 

Given a statistic stat on ©„, the distribution of stat over ©„ is the polynomial 

When a = 1", we will simply write Two statistics, say stat on Obj and 

stat' on Obj', are said to be equidistributed if their distributions are equal. One 
particularly nice way to prove equidistribution is to give a bijection / : Obj —>■ Obj' 
such that stat'(/(cr)) = stat(cr) for every a G Obj. Our main result will be a 
bijection of this form. Finally, it will be convenient to use interval notation for the 
integers, i.e. [a, b] is the set of integers at least a and at most b. 

In [Macl5], MacMahon showed that inversion number and major index are 
equidistributed over &a, and that 

DTiq) = DTKq) = 

where we use the standard g-analogs 

[n]q\ = [n]q[n-V\q...[\]q [n]q = . 

MacMahon’s proof was not bijective; the first bijective proof of this fact was given 
in [Foa68]. A second proof, essentially due to Carlitz [Car75], is sometimes known 
as the insertion method. It will be the template for our bijections in Section 3. 
Part of our main result^ can be viewed as a bijective proof of the identity 

des((7) 

(I) ^ n (1+= E n (i+^/g”"°''+') 

i—1 jGDes{cr) 

for any composition a. The case a = I” was proved in [RWI5]. In order to provide 
our bijective proof, we will rephrase this statement as an equidistribution result on 

^To state our main result in full we will need more notation, which we will define later. 


Gll, 0^25 ■ ■ • 5 q ■ [ 0 ^ 2 ]^ ■ ■ • • [o^njg- 
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ordered multiset partitions. Then we will generalize Carlitz’s insertion method to 
ordered multsiet partitions in order to complete our proof. 

1.2. The insertion method for 6„. One consequence of MacMahon’s equidistri- 
bution theorem is a pair of recursions for the distributions of the inversion number 
and the major index over ©„: 

( 2 ) DT{q) = [n],D^lM DTKq) = [n],D^^M- 

On the other hand, these two statements imply MacMahon’s result. Carlitz’s in¬ 
sertion method gives bijective proofs of these statements which can be combined 
to build a recursive bijection -ipn '■ &n &n such that maj(i/’(cr)) = inv(CT). We 
say that ipn maps the inversion number to the major index. We outline Carlitz’s 
insertion method below. 

To prove the left statement in (2), one simply considers all the possible ways 
to insert n into a permutation in &n-i to create a permutation in 6n- It is clear 
that, for cr G &n-i, inserting n after the first i elements of a creates n — i — \ new 
inversions and does not affect the previously existing inversions. For example, for 
<j = 5167324 G © 7 , we can “label” these positions with subscripts that give the 
number of new inversions created by inserting an 8 at that position: 

756156473822140 . 

This proves the inversion side of (2). The key to the insertion method is that 
something similar is true for the major index. In particular, we can label the 
spaces between elements of ct G ©n-i, along with the left and right ends, according 
to the following scheme: 

(1) Label the position after cr„_i with a zero. 

(2) Label the descents of a right to left with 1,2,..., des(CT). 

(3) Label the position before u\ with des(CT) -I- 1. 

(4) Label the ascents of a from left to right with des((T) -I- 2,..., n — 1. 

For example, cr = 5167324 receives the following labels in this setting: 

453156672812740 

These labels give the change in major index that comes from inserting n at that 
position; one proof of this fact can be found in [HLR05]. This completes the proof 
of (2) and also gives a bijection that takes the inversion number to the major 
index. We include an example of this bijection in Figure 1. To compte i/’5(52143), 
we remove the 5 and count the number of inversions lost by removing 5. In this 
case, we have lost 4 inversions. We record this number in the third column and 
the resulting permutation in the a column. We repeat this process until we have 
reached n = 1 and filled the first three columns of the table. To build our new 
permutation, we recursively place n at the position that receives label i in the 
major index labeling. These labels have been italicized in the example. 


1.3. The insertion method on ©„. It is natural to hope that this proof can be 
extended to permutations that may contain multiple copies of the same number. 
That is, we would like to give insertion proofs that 


( 3 ) 


DTiq) 





DT'^iq) 



DT-\q). 
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n 

a 

Change in inv 

V’n(cr) 

5 

52143 


24153 

4 

2143 

4 

2234i1^3o 

3 

213 

1 

22jl33o 

2 

21 

0 

22i1o 

1 

1 

1 

jIq 


Figure 1. We compute V'5(52143). 


where a~ = (oi, 02 , • ■ • > an-i)- Such proofs would imply MacMahon’s equidis- 
tribution theorem and provide a bijection between inversion number and major 
index. 

The inversion side cooperates nicely. As before, inserting an n to the right of i 
elements of cr G ©„ increases the inversion number by \a~ \ —i. Hence this insertion 
can create between 0 and |q;“| inversions. Furthermore, the position of a new n 
has no affect on the number of inversions added by other n’s; in other words, each 
insertion is independent of the other insertions. This allows us to compute {q) 
from 


DT{q)=D^l{q) I n 


2=0 


1 — q^x 







jaf 

1 









q 


To prove the major index side of (3), we essentially recreate the bijection con¬ 
structed in [FH08, CPYYIO]. Let (j^f j) denote the family of fc-element multisets 
containing elements drawn from the set S. We would like to establish a bijection 


X 


[0,|a-|] 


6c 


such that 


maj (C‘'^(cr, B)) = maj(cr) + ^ 

b&B 

Such a map would provide a combinatorial proof of the major index side of (3). 
Before inserting any n’s, we label a G &a- in a manner reminiscent of Section 1.2: 

(1) Label the position after cr|c-| with a zero. 

(2) Label the descents of a right to left with 1,2,..., des(CT). 

(3) Label the position before u\ with des(CT) -I- 1. 

(4) Label the non-descents of a from left to right with des(CT) -|- 2,..., |q;“|. 

Write B = {bi > b 2 > ■ ■ ■ > We insert an n into the position labeled &i. 

Then we go through the labeling process again, stopping once we have used the 
label bi- We insert an n into the position labeled 62 - We repeat this process until 
we have processed each element of B. We omit the proof that this map satisfies the 
desired properties, which can be found in [FH08, CPYYIO]. Instead, we will work 
through an example. 

Let a = {2,1, 3, 2}, cr = 323113 G 6a-, and B = {5^}. We note that maj(CT) = 4. 
We begin by labeling cr according to the labeling associated with the major index. 
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a 

a 

B 

i^Aicr) 

11,2,1,3} 

2443214 


4432124 

4 4382211 2(7 
4382211 2 o 

{ 1 , 2 , 1 } 

2321 

{3,3,0} 

3322il42o 

{ 1 , 2 } 

212 

{ 2 } 

221 I 320 

2 ilo 

{ 1 } 

1 

{ 1 , 0 } 

jIq 


Figure 2. An example of the map 4’a for a = {1, 2,1,3}. 

We place a 4 at the label 5 to get 3231413. We relabel this permutation, stopping 
when we use the label 5. 


4332532 I 4il 3o 

Then we insert a 4 at the position labeled 5 to get 32431413. As desired, 

maj(32431413) = 14 = maj((T) + ^ 5 = 4 + 5 + 5. 

6eB 

Just as before, these insertion maps can be combined to yield a bijection ipa '■ 
&a ®a that takes inversion number to major index. We illustrate tpa with the 
example in Figure 2. As in Section 1.2, we fill the first three columns of the table 
from top to bottom by removing all copies of the largest element and recording the 
multiset of inversions lost during each removal, which we call B. Then we fill the 
fourth column by using the labeling associated with the major index to repeatedly 
insert a new element at the position that received the largest remaining label in B. 

2. Ordered Set and Multiset Partitions 

2.1. Definitions. The ordered set partitions of order n with k blocks are partitions 
of the set {1,2,..., n} into k subsets (called blocks) with some order on the blocks. 
We write this set as OVn,k- For example, 13|45|2 G 07^5,3, where we have listed 
each block as an increasing sequence and we have used bars to separate blocks. It 
is not difficult to see that OVn,n = &n, so ordered set partitions are a natural 
extension of permutations. 

More generally, given a composition a of length n, the ordered multiset partitions 
OVa,k are the partitions of the multiset A{a) = : 1 < z < n} into k ordered 

sets, which we still call blocks. For example, 24|134|2 G OV(^ 1 ^ 2 , 1 , 2 ), 3 - Note that, 
although we are dealing with the elements of a multiset, each block is still a set. 
The analogous objects where blocks are also multisets will not arise in our work. 

So far, we have written each block of an ordered set or multiset partition in 
increasing order from left to right. We will often wish to use the opposite notation, 
i.e. we will write each block in decreasing order from left to right. Furthermore, we 
will use stars as subscripts to “connect” elements in the same block instead of bars 
to separate blocks. For example, the ordered multiset partition 24|134|2 is written 
as 4*2 4*3*1 2 in this new notation. We will refer to an ordered multiset partition 
written this way as a descent-starred permutation, since every permutation of the 
given multiset with some (but maybe not all) of its descents “starred” corresponds 
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to an ordered multiset permutation in this fashion. More formally, we define the 
descent-starred permutations of A(oi) with k stars as follows: 

K,k = S):aG6A,SC Des(a), |5| = k}. 

The set S corresponds to the entries of cr which are followed by stars. Then there is 
a straightforward bijection OVa,k ^ given an ordered multiset partition, 

we write its blocks in decreasing order from left to right, add stars between adjacent 
elements that share a block, and remove the bars. 

2.2. Statistics. We will study the following four statistics on ordered multiset 
partitions. All four of them appear in connection with a certain operator in the 
theory of diagonal harmonics which we call the Garsia-Haiman delta operator. 
Two of them (inv and maj) are directly involved with the statement (1) given in 
the introduction. 

First, given tt G 0'Pa,k, inv(7r) counts the number of pairs a > b such that a’s 
block is strictly to the left of 6’s block in tt and b is minimal in its block in tt. We 
call these pairs inversions. For example, 15|23|4 has two inversions, between the 5 
and the 2 and the 5 and the 4. 

For any tt G OVa,k, we number tt’s blocks tti, 712 ,..., from left to right. Let 
TT^ by the hth smallest element in tt^, beginning at h = 1. Then the diagonal 
inversions of tt, written Dinv(7r), are the triples 

{{h,i,j) --I <i < j < k, 7rf > tt'*} U {{h,i,j) ■ ^ <i < j < k, 7rf < 7r''+^}. 

The triples of the first type are primary diagonal inversions, and the triples of the 
second type are secondary diagonal inversions. We set dinv(7r) to be the cardinality 
of Dinv(7r). For example, consider the ordered multiset permutation 24|134|2. It is 
helpful to “stack” the elements in each block vertically, obtaining the diagram 





4 


2 


to 


The primary diagonal inversions are (1,1,2) (between the leftmost 2 and the 1 in 
the first row) and (2,1, 2) (between the 4 and the 3 in the second row) and the only 
secondary diagonal inversion is (1,1, 2) (between the leftmost 2 in the first row and 
the 3 in the second row), for a total of three diagonal inversions. 

To define the major index of tt, we consider the permutation cr = cr(7r) obtained 
by writing each block of tt in decreasing order. Then we recursively form a word 
w by setting wq = 0 and Wi = Wi-i + is minimal in its block in tt) for each 
i > 0. Then we set 

niaj(7r) = ^ Wi. 

i: <Ti>CTi+i 

Using the ordered multiset permutation tt = 24|134|2 again, we obtain a = 424312 
and w = 0011123, beginning with wq = 0. The descents of a occur at positions 1,3, 
and 4, so maj(7r) = iwi + rua + r (;4 = 0 + 1 + 1 = 2. 

There is an alternate definition of the major index which we will use in some 
of our proofs. It is clear from the definition above that for any tt G OVn,m the 
major index defined here is equivalent to the major index defined on permutations 
in Section 1. Now we consider what happens to maj(7r) if we decide to insert 
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a star after a descent at position d. We note that Wi decreases by 1 for each 
i > d. Therefore the major index of tt has decreased by 1 for each descent weakly 
to the right of position d. Therefore, if {a, S) is the descent-starred permutation 
representation of tt, we can write 

(4) maj(7r) = maj(CT) — ^ | Des(CT) n {i, i -|- 1,.. .}|. 

ies 


This shows that 

des((T) 

n n (i+w) 

*=1 z\a\-h 

Finally, we define the minimum major index of tt as follows. We begin by 
writing the elements of tt^, in increasing order from left to right. Then, recursively 
for 7 = fc — 1 to 1, we choose r to be the largest element in tt^ that is less than or 
equal to the leftmost element in Tr^+i, as previously recorded. If there is no such r, 
we write tt^ in increasing order. If there is such an r, beginning with tt^ in increasing 
order, we cycle its elements until r is the rightmost element in TTj. We write down tt^ 
in this order. We continue this process until we have processed each block of tt. For 
example, consider the ordered multiset permutation tt = I3|23|I4|234. Processing 
the blocks of tt from right to left, we obtain 3I234I234. We consider the result as 
a permutation, which we denote t = T(7r), and define 


minima j(7r) = i 

i: Ti>Ti+i 


i.e. the major index of the permutation t. The name minimaj comes from the fact 
that minimaj (tt) is equal to the minimum major index achieved by any permutation 
that can be obtained by permuting elements within the blocks of tt. 


3. EQUIDISTRIBUTION ON ORDERED MULTISET PARTITIONS 


In this section, we prove the following equidistribution theorem for ordered mul¬ 
tiset partitions. 


Theorem 1. For any composition a, 


The reader may have noticed that ( 9 ) is not included in the list of 

equidistributed polynomials; that is because we have not proved this case at this 
point. We describe the unique difficulties of this case in Subsection 3.4. 

Our proof of Theorem 1 is bijective and employs a generalization of Carlitz’s 
insertion method from permutations to ordered multiset partitions. We describe 
“insertion maps” for inv, dinv, and maj in Subsections 3.1, 3.2, and 3.3, respectively. 
These maps will be of the form 


Cm ^ X 

Cm : X 


[0,£-I] 

an — k + £ 

[0,£-l] 

an — k + £ 



X 

X 

X 



OVa,k 

OVa,k 

OVa,k 
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where we have used some new notation here. For any composition a of length n, 
a~ is the composition obtained by removing the rightmost entry in a. For a set 
S', (f) is defined to be the collection of all subsets of S of size k. Similarly, (j^f |) 
is the collection of all multisets of size k whose elements are taken (possibly more 
than once) from S. By “insertion maps,” we mean that they satisfy the properties 

inv ^)) = inv(7r) + '^u+'^b 

ueu beB 

U,B)^ = maj(7r) + ^ u + ^ 5 
uGU beB 

dinv ^)) = dmv(7r) + “ + XI 

ueu beB 

When k = \a\ these maps will reduce to the insertion processes defined in Subsection 
1.3. We will use these maps to construct a bijective proof of Theorem 1. Subsection 
3.5 contains more information about the shared distribution of the polynomials in 
Theorem 1. 

Finally, it is natural to wonder if our results can be transferred to ordered parti¬ 
tions of a multiset into multisets instead of sets. For example, 113|23|1 is one such 
object. At this point, we have not managed to accomplish this task; in particular, 
the major index statistics seems to behave differently in this setting. However, 


3.1. Insertion for inv. Recall that we need to define a map of the form 


[0,f-l] 




[0,f] 


or 




We can think of the set U G providing the increases in inv that come 

from adding a new n without creating a new block, and the multiset B G 
as providing the increases in the statistic that come from adding a new n while 
creating a new block. These maps will cooperate with our inversion statistic in the 
following manner: 


(5) inv AAkA'^^ ^)) = inv(7r) + X “ + X ^ 

ueU bGB 


Given tt € we label each block of tt from right to left with the numbers 

0,1,..., f — 1. We repeatedly remove the largest element from the multiset U U B, 
taking the element from U if the largest elements are equal. Call this element i. 
If i came from U, we place an n in the block that received the label i. If i came 
from B and is equal to £, we place an n as a new block to the left of the block 
that received the label i. Hi came from B and is less than £, we place an n as a 
new block just to the right of the block labeled i. The resulting ordered multiset 
partition is (tt, U, B). 

For example, say a = {I, 2, 2,3}, fc = 6, and f = 5. We consider tt = 3|1|2|2|I3 S 
or a-, 5 , U = {2,0}, and B = {3}. Then U VJ B = {3,2,0}. We label tt as follows, 
with the labels written as subscripts at the end of each block: 


34|l3|22|2i|I3o 
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The largest element in [/ U i? is 3 and it comes from i?, so we insert a 4 at the 
position labeled 3 as a new block to the right of the block labeled with the 3. 

34|l3|4|22|2i|13o 

Now the largest remaining element of 17 U B is 2 and it comes from 17, so we put a 
4 into the block labeled 2. 

34|l3|4|242|2i|13o 

Finally, we insert a 4 into the block labeled 0 to obtain (tt, U, B). 

34|l3|4|242|2i|134o 
We can check that (5) holds here. 

10 = inv(3|l|4|42|2|431) 

= inv(3|l|2|2|31) + 

«e(7 bGB 

= 5 +(2 + 0)+ 3. 

Lemma 3.1.1 (Insertion for inv). For any composition a of length n and positive 
integers £ < k, (f'ffl. ^ is well-defined and injective. The image of ^ is the 
ordered multiset partitions it G OVa,k with exactly k — i singleton blocks containing 

n. Furthermore, for any tt G OVa,i, U G and B G 

inv (C,fc, U, B)) = inv(7r) + E “ + E 

ueu bes 

Proof. The only way ^ could not be well-defined is if we tried to insert two n’s 
into the same block. Since 17 is a set, this does not occur. The statement about 
the image of i follows from the definition. Furthermore, the insertion map is 
clearly bijective, and any function that bijects onto its image is an injection. 

We will prove that whenever we remove the largest element i from U U B and 
introduce a new n to our ordered multiset partition as described in the insertion 
map, we introduce i new inversions to the ordered multiset partition. Specifically, 
we create inversions between this new n and the minimal elements of the i labeled 
blocks of TT that are to the right of the block that received label i. Since n is 
the largest entry, we do not create any new inversions that end at n. Finally, we 
note that we do not destroy any inversions that existed before we inserted this new 
n. □ 

3.2. Insertion for maj. To define we will view tt G as a descent- 

starred multiset permutation (cr, S). We will label the unstarred positions of a 
as in Section 1. Specifically, we label the unstarred descents from right to left, 
then the non-descents from left to right, using increasing labels 0, !,...,£. Let 
t/+ = {m + 1 : M G 17}. We repeatedly remove the largest element i from U~^ U B, 
taking i from B if the largest elements are equal. Then we proceed through the 
following algorithm: 

(1) Insert an n at the position labeled i. 

(2) Move each star to the right of the new n one descent to the left. 

(3) If i came from [/+, star the rightmost descent. 
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(4) Relabel as before, stopping at the label i i came from B and i — 1 \i i 
came from U'^. 

When we have used each element of 17+ U B, the result is S), U, B). 

For example, let us again consider a = {1,2,2,3}, k = 6, and £ = 5 with 
(ct,S') = 3 1 2 2 3*1 G OVa-, 5 , U = {2,0}, and B = {3}. Then C/+U R = {3,3,1}. 
We label (cr, S) as follows. 

23il324253*lo. 

We insert a 4 at the position labeled 3 and then move all stars to the right of that 
position one spot to their left. Since we took 3 from B, we do not star the rightmost 
descent, resulting in 3 1 4*2 2 3 1. We re-label to obtain the following. 

332 I 4*2 2 3ilo 

Again we choose the position labeled 3. This time we star the rightmost descent 
after shifting stars because 3 came from 17+, yielding 4 3*1 4 2 2 3*1. Finally we 
label this element 

423*1 4i2 2 3*lo. 

We insert a 4 at the position labeled 1, shift stars, and star the rightmost descent 
to get 4 3*1 4 4*2 2 3*1. We check that this new permutation has the desired major 
index. 

10 = maj (4 3*1 4 4*2 2 3*1) = l + l + 3 + 5 
= maj (3 1 2 2 3*1) 4- 

ueU bGB 

= {l + A) + {2 + 0) + 3. 

Now we prove that this process satisfies the necessary properties. 

Lemma 3.2.1 (Insertion Lemma for maj). For any composition a of length n and 
positive integers I < k, ^ is well-defined and injective. Furthermore, for any 

n G U G and B G ((L°i)), 

maj u, B)^ = maj(7r) -t ^ M -f ^ b. 

ueu beB 

Proof. To show that (j\Ak i is injective, we describe its inverse. Given some descent- 
starred multiset permutation (t, T), we first check if the rightmost descent of (r, T) 
is starred. If it is, we remove the star and prepare to add an element to U. Other¬ 
wise, we prepare to add an element to B. We scan t for the rightmost n which is 
either at the right end of r or between two entries such that the entry to the left 
of n is greater than the entry to the right of n. If there is no such n, we choose the 
leftmost n in r. We move all stars that are weakly to the right of this n’s position 
one descent to their right and then remove n. Say that, at this point, we have 
decreased the original major index of (t, T) by i. We add i to either U or B, as 
decided above. Then we repeat this process until we have removed all n’s. It is 
important to note that the inverse does not depend on knowledge of £; therefore, 
each p G OVa,k is in the image of e. ^ unique value of 1. 

In the remainder of the proof, we show that the map cooperates with the statistic 
maj as proposed in the lemma. We consider the labeling of the descent-starred 
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multiset permutation (cr, S) equivalent to tt at any step during the insertion process. 
Say i is currently the largest element of [/+ U B. 

Assume first that the position labeled i is a descent. We use d to denote the 
number of starred descents to the right of this position. By the insertion method 
discussed in Subsection 1.3, inserting n into the position labeled i increases maj(cr) 
by i+d. Furthermore, we have not created a new descent, so the number of descents 
weakly to the right of any starred position has remained the same. Therefore, by 
the alternate definition of maj(7r) in (4), we have increased maj(7r) hy i + d after 
Step 1. 

For Step 2, we move all stars to the right of the position labeled i one descent 
to their left. Since position i contains an unstarred descent, this is always possible. 
Furthermore, each of these d stars have picked up an additional descent that is 
weakly to their right. Using (4) again, we see that the change in maj(7r) after Step 
2 is i + d — d = i. 

Finally, we need to consider if i came from U~^ or B. Hi came from C/'*", we 
star the rightmost descent. This subtracts 1 from maj(7r). In either case, we have 
increased the major index of tt by the amount equal to the element from U or B 
corresponding to i. 

By the insertion process from Subsection 1.3, we can relabel the resulting descent- 
starred permutation and repeat the process as long as we bound the labels as 
described in Step 4. Then, by the same argument as we used above, the insertion 
process will modify the major index as described in the statement of the lemma. 

Now we consider where the argument must change when the position labeled i is 
not a descent. We still use d to denote the number of starred descents to the right 
of position i, and we set c to be the number of starred descent to this position’s left. 
Since every starred descent occurs before the position labeled i in the labeling order 
for maj(cr). Step 1 increases maj(7r) hy i + c + d. For Step 2, the position labeled 
i still contains an unstarred descent, so we can still move the stars as described. 
As before, this means that each of the d stars to the right of the position labeled i 
adds a descent to its right, contributing —d to maj(7r). Furthermore, inserting n at 
a non-descent creates a new descent, so each of the c stars to the left of the position 
labeled i has added a descent to its right, contributing —c to maj(7r). Therefore 
the total increase of maj(7r) is i + c + d — c — d = i. Steps 3 and 4 do not depend 
on whether we are inserting at a descent or a non-descent. □ 


These two insertion maps work together to provide a bijection 'tpa,k ■ OVa,k —> 
OVa,k that takes inversion number to major index. The bijection is described 
recursively as follows. 

(1) Given an ordered multiset partition p G OVa,k, choose ^ such that p has 
k — £ singleton blocks containing n. 

(2) Set (tt, U, B) to be the inverse of p under 

(3) Recursively send tt to tt' = 

(4) Set ■i/’a,fc(p) = U, B). 


Finally, in order to begin the recursion, we declare that tpim jyi is the identity map. 
We work through an example of this bijection in Figure 3. 
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a 

k 

IT 

U 

B 

V’A.fc(7r) 

{3,1,2,!} 

3 

134|1|3|12 



3*2*1 14 3 1 

{3,1,2} 

2 

13|1|3|12 

{4} 

0 

23*2*l3l43ilo 






i2*l2l 3*lo 

{3,1} 

1 

1|1|12 

{2} 

{1} 

l2*l2l3lo 

{3} 

0 

1|1|1 

{0} 

0 

J 121310 


Figure 3. An example of the map 'tpA,k{T^)- 


Proposition 3.2.1. For any composition a and positive integer k, '4’a,k is a bijec- 
tion with the property 

niaj(i/>a.fe(/i)) = inv(p) 

for any p € OVa,k- 

Proof. We will work by induction on n, the length of a. If n = 1, Then the multiset 
is {1“^} and k must be equal to oi. In this case OVa.k only has I element, which 
has oi parts all equal to 1. This element clearly has inv = maj = 0. We defined 
ipa,k so that it is the identity in this case, which clearly is a bijection and satisfies 
maj(i/’Q,fc(7r)) = inv(7r) for the unique tt G OVa,k. 

If n > 1, take any element p € OVa.k- We choose t, tt, tt', U, and B as instructed 
in the definition of ipa.k- The images of ^ for £ = I to A: partition OVa,k into 
the subsets consisting of elements which have fc — t singleton n blocks. (Since we 
assume each ai > 0, an element cannot consist entirely of singleton n blocks.) We 
also noted while proving Lemma 3.2.1 that each p € OVa.k is in the image of (f^’^k i 
for a unique value of Since each of these insertion maps is invertible, i/'a.fc is a 
bijection. 

Finally, we use Lemmas 3.1.1 and 3.2.1 along with the inductive hypothesis to 
compute 

niaj('i/'a.fe(p)) = maj(7r') + '^u + '^b 

ueU bGB 

= inv(7r) + u + b 
ueu bGB 

= inv(p). □ 

We work through an example of the map 'ipA.ki'^) in Figure 3. We repeatedly 
remove all of the largest elements (and their stars) from the starred permutation 
and recording the number of inversions lost in the U and B columns. Starred 
elements contribute to the U column and unstarred elements contribute to the B 
column. Once we have reached the final row, we use this information to build the 
'f’a,k{{(k, S)) column from bottom to top. We use the elements of C/+ U M to select 
the positions at which to insert new largest elements. This insertion follows the 
procedure laid out in the definition of 

There are a number of other consequences of our proof. For example, the right- 
to-left minima of a permutation cr are the entries at such that, for all j > i, 
Ui < Uj. In [RW15], Remmel and the author proved that the a = I" case of 
'f’a.k preserves the right-to-left minima of ct, where (cr, S) is the descent-starred 
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permutation representation of an ordered partition of a. The same is true for 
general ipa^k- 

Corollary 3.2.1. For any element of tt G OVa,k considered as a descent-starred 
permutation consider its image 4’a,ki'^) as the descent-starred permutation 

(r, T). Then a and r have the same right-to-left minima. In particular, an = Tn- 

The crux of the proof is that the insertion algorithms only change the last element 
of a when 0 is an element of the multiset B. 


3.3. Insertion for dinv. In order to prove that dinv is equidistributed with inv 
and maj, we define an insertion map for dinv 


/'o.M : X 


a 


[ 0 ,^- 1 ] 

-k + I 


[ 0 ,£] 

k-I 


or. 


,k' 


Given tt € Ora-,e, we will use two different labelings to insert the n’s into tt. We 
label the positions between the blocks, as well as the positions at either end of tt, 
with the labels 0,1, from right to left. We will call these the gap labels. 

We will label the £ blocks of tt with the labels 0,1,... ,£ — 1. Set h to be the 
maximum height of any element in tt. We obtain the block labels of tt by, beginning 
with the label 0, repeatedly 

(1) labeling the elements of tt at height h from left to right with increasing 
labels, and 

(2) decrementing h. 

We repeat until h < 0. For example, if tt = 124|2|13|134|1, the block labels of tt are 
0|3|2|1|4. 

With these labels in hand, we define O, B) by inserting an n into each 

block that receives a block label u € U and into each gap that receives a gap label 
b G B. As usual, the key is to prove that this map cooperates with the statistic 
dinv. 


Lemma 3.3.1 (Insertion for dinv). For any composition a of length n and positive 
integers £ < k, is well-defined and injective. The image of is the 

ordered multiset partitions tt G OVa,k with exactly k — £ singleton blocks containing 

n. Furthermore, for any tt G OVa,i, U G and B G 

dinv U, B)) = dinv(7r) + ^u-G^b. 

uGU bGB 

Proof. It is clear that inserting an n at a gap labeled b creates b new diagonal 
inversions, one with each block to the right of the gap, and does not affect any 
other diagonal inversions. 

Now say we insert an n into a block labeled u. Say that tt^ had size s before 
we added an n. After inserting an n, it has size s + I with an n at height s + I. We 
claim that we have created one new diagonal inversion for each block of size greater 
than s and for each block of size s that is to the left of TTj. First, consider a block 
TTj with \TTj\ > s. If j > i, then (s + l,i,j) is a new primary diagonal inversion; 
if j < i, then (s,j,i) is a new secondary diagonal inversion. There are no other 
new diagonal inversions between and ttj. Now we consider ttj with j < i and 
\TTj\ = s. There is a new secondary diagonal inversion (s,j,i). By the definition 
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of our insertion map, there are exactly u such blocks, so we have created u new 
diagonal inversions. 

Finally, we note that, since two n’s cannot form a diagonal inversion, each in¬ 
sertion of an n does not affect previous or subsequent insertions. □ 

It follows from Lemma 3.3.1 that We can form 

a bijection between any pair of these statistics by using the definition of ipa,k as a 
template. Furthermore, any such bijection preserves right-to-left minima. 

3.4. The statistic minimaj. Data computed in Sage suggests that minimaj shares 
the distribution of inv, maj, and dinv; unfortunately, we cannot prove this with the 
techniques currently available to us. However, data yields the following conjectures 
as to how this statistic relates to the previous statistics. Given an ordered set 
partition tt G OVn, we say the shape of tt, written shape(7r), is the composition 
whose zth block is equal to the size of tt*, the zth block in tt from left to right. 

Conjecture 1. For any composition (3 of length n, 

Y. 9'"”"' = E ^minimaj(7r) 

T^eOVn T^eOVn 

shape(7r)=/3 shape(7r)=/? 

Also, for any composition a 

D~^'\q) = D^liq) = D^^{q) = Dtj{q). 

Note that neither statement directly implies the other. 

In personal communication, Brendon Rhoades has notified the author that he is 
currently preparing a proof of both parts of this conjecture. 


3.5. The Mahonian distribution on 0'Pa,k- In this subsection, we describe 
the distribution shared by the statistics inv, maj, and dinv on OVa.k- Define the 
Mahonian distribution on OVa.k to be the polynomial 

Do.,k{q) = D“X(g) = D^^{q) = Ddj(g). 

We know from MacMahon’s theorem that Dq,_|q|(( 7 ) = ] . In general, we 

can only give a recursive description of Dct,k{(l)- Applying standard q-binomial 
identities to the insertion maps given above, we see 


( 6 ) 


-htct,k(,Q) 





£ 


'k 

2 J 

an — k + £ 

q 

_£_ 


A^) 


with initial condition 


(7) 7)(ai),fc(9) = x(fc = ai)- 

We can cancel terms of (6) to obtain the identity 

k 

( 8 ) Da,k{q) = ^ ^ 

/=1 


a-n — k + k — a-n^ k — i 


Da-Aq)- 


We can obtain another expression for this polynomial in the special case ai = 
.. = «„ = a. Before we can state this expression, we must define a g-analog of the 



EQUIDISTRIBUTION ON ORDERED MULTISET PARTITIONS 


15 


(generalized) Stirling numbers of the second kind. The (? = 1 case of these poly¬ 
nomials appear in [BPS03], Equations (20) and (21). We define these polynomials 
recursively by 


(9) 

( 10 ) 


= Xik = a). 


a — k + i 








-4! 


Note that, at a = 1, the recursion simplifies to the g-Stirling numbers Sn,kiq)- 


Proposition 3.5.1. When a = a", 

When 0 = 1, this formula reduces to the formula 

h^n,k{.q) — [k]q^-Snjk{q) ’ 

Proof. We work by induction on n. When n = 1, the right-hand side of Proposition 
3.5.1 equals 


= “> 

which is equal to the left-hand side by (7). 

If n > 1, we use the induction hypothesis with the recursion (8) to compute 


D. 






k 

o — fc-l-£, fc — o, fc — . 


(H,) 




[k]ql 


r__ q(“) roi 

[a-k + e]ql[k-a]ql[k-e]qr^-^^^^^^ 


(H,!) 

[k]q'- 






a-k + e\\k-eA 


^ {k]q'- c(a)(- \ 

([o],!)” 


by (9) with i = £. 


□ 


It would be interesting to give a more combinatorial proof of Proposition 3.5.1, 
especially one that would shed light on why each of the three terms [fc]g!, ([a]g!)”, 
and Skj^iq) appears. 


4. Extending Macdonald polynomials 

In this subsection, we apply our inv and maj statistics to the combinatorial 
definition of Macdonald polynomials for hook shapes, as given (for any shape) in 
[HHL05]. This yields functions whose coefficients are four-variable polynomials 
instead of the usual two-variable polynomial coefficients. Our previous work allows 
us to prove that these polynomials are symmetric and to expand them into Schur 
functions. 
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For convenience, for any statistic stat, let 


stat[a,6] (cr) = Stat(craO'a+l • • • df,). 

We will also need two new statistics, one on permutations and one on descent- 
starred permutations: 

rlmaj(CT) = (n — i) 

iG'Des{(7) 

rlmaj {{a, S)) = rlmaj( ct) - ^ | Des( ct) n [1, i]|. 

ies 

We set Hn,mix; q,t,u,v) equal to 


E 

cre{l,2,...} 




X n 

zGDes((7)n[m+l,n] 


1 + 


1 (^) + l 


I Des((7)n[l,m] I 

n a 

i=i 


,/t^) . 


We will refer to these polynomials as starred Macdonald polynomials of hook shape. 
These polynomials can be thought of as a sum over all descent-starred multiset 
permutations {a, S) in the Young diagram for the shape (n — m, 1™) where we cal¬ 
culate our maj and inv statistics down the column and across the row, respectively. 
When u = V = 0, we obtain the (modified) Macdonald polynomial for the shape 
(n — m, 1™), as proven in [HHL05]. Here is an example filling for n = 8 and m = 3. 




6 * 

m 

m 

4 

00 


The numbers followed by stars are the starred descents. The weight of this Filing 
would be q'^t, since the descent-starred permutation 6*1748 has 2 inversions, both 
ending at the 4, and 25*36 has major index equal to 1. 

Our main result in this section allows us to transfer many important properties 
of the Macdonald polynomials to the starred Macdonald polynomials of hook shape. 
The u = V = 0 case of this result was originally proved in [Ste94]. The result requires 
some definitions on standard Young tableaux. The descent set of a standard Young 
tableau T (in French notation) is the set of all i such that t -I-1 is strictly north (and 
weakly west) of i in T. Then, for any standard Young tableau T with n entries. 


maj(T) = i 

ieDes(T) 


rlmaj (T) = (n — i). 

'iGDes(T) 


For example, here is a standard Young tableau with descent set {2, 5}. 


3 

4 

6 


1 ” 

to 


7 


Theorem 2. The starred Macdonald polynomials of hook shape are symmetric. 
Furthermore, for A F n the coefficient of the Schur function s\ (x) in 
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Hn,m{x;q,t,u,v) is equal to 
TeSYT(A) 

I Des(T)n[m+l,n] I | Des(T)n[l,m] | 

X (l + uq~^) (1 + vt~^) . 

i=i j=i 

Proof. We begin with a descent-starred multiset permutation {a, S). Then we apply 
the bijection 

7 = complement o reverse o o reverse o complement 

to the descent-starred multiset permutation {am+i ■ ■ ■ cTn, S r\[m +1, n]) for suitable 
l3,£. Since preserves the rightmost letter, 7 preserves ak+i and sends the inv of 
(um+i ■ ■ ■ CTn, S n [m + l,n\) to the rlmaj of the resulting descent-starred multiset 
permutation. Hence, equals 

(11) V ^<T^rImaj[j,^i_„](<T)^maj[i_j,] (o' 


( 12 ) 


E 

I Des((T)n[/i:+l,n] I | Des((7)n[l,fc] | 

X (l + u/q^) J]^ (l + v/P). 

i=i j=i 


For any composition a, the coefficient of in this expression is just the sum over a 
that are permutations of the multiset {1“^, 2“^ ,•••}• We would like to show that the 
coefficients of x°‘ and a;“ are equal, where is obtained from a by switching 
and ar+i- To do this, we apply a procedure known as r-pairing to the permutation 
T. We illustrate r-pairing via the example in Figure 4. To perform r-pairing on 
a sequence, we begin by temporarily ignoring all entries not equal to r or r -f 1. 
Then we pair off adjacent occurrences of r -f 1 and r, ignoring previously paired 
entries and iterating this pairing procedure. When there are no more such pairs, 
we replace each un-paired occurrence of r with an r -I- 1 and vice versa. Finally, 
we re-insert the entries we had temporarily removed in their initial positions. We 
provide an example above with r = 2. 

For our purposes, it is enough to know that r-pairing replaces cr with a permuta¬ 
tion of the multiset {l“b ..., r“’’+'^, (r -I- 1)“’",..., }, and that this permutation has 
the same descent set as a. Therefore r-pairing does not alter any of the expressions 
in (11), so Hn^m is symmetric. 

Furthermore, if we apply the Robinson-Schensted-Knuth correspondence [Sta99] 
to each permutation cr involved in the coefficient of x^ for a partition A, we see that 
the coefficient of x^ in Hn^m is 

rlmaj[,„+i „](T) „](T) 


K 


q 


A,(n—m,l^) / ^ 

TgSYT(A) 

I Des(T)n[m+l,n] I | Des(T)n[l,m] | 

X Y[ (l + u/q^) (l + v/P). 

i=l j=l 

Here is the Kostka number. Translating from the coefficient of x^ to 

the coefficient of SA(a;) exactly consists of removing this Kostka number, so the 
theorem follows. □ 
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24231243331324123321 
2 23_^ 333 ^ 23^ 

2 2 333_23 

2 2 33 3 

3 3 22 2 

3 3 223_22 

3 31^ 223 ^ 22^ 
34331242231324122321 


Figure 4. An example of r-pairing with r = 2. 


One consequence of Theorem 2 is the identity 

(13) Q: 1^) — 1^5 R)- 

This parallels the well-known identity 

(14) Hf,{x;q,t) = Hfj,^{x;t,q) 

for Macdonald polynomials. For Macdonald polynomials of hook shape, there is a 
direct bijective proof of (14). We would like to find a similar bijective proof of (13). 

In future work, we hope to define and explored starred Macdonald polynomials 
for non-hook shapes. The full combinatorial formulation of Macdonald polynomials 
in [HHL05] essentially allows one to either star horizontally or vertically, giving 
analogs of Macdonald polynomials with three-variable polynomials for coefhcients. 
At this point, we have been unable to define a four-variable generalization in the 
non-hook case that retains desirable properties such as symmetry. It would also be 
interesting to develop connections between these polynomials and Garsia-Haiman 
modules. 
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